Loop cosmological implications of a non-minimally coupled scalar field. 
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Non-minimal actions with matter represented by a scalar field coupled to gravity are considered 
in the context of a homogeneous and isotropic universe. The coupling is of the form —^^(f>^R. 
The possibility of successful inflation is investigated taking into account features of loop cosmology. 
For that end a conformal transformation is performed. That brings the theory into the standard 
minimally coupled form (Einstein frame) with some effective field and its potential. Both analytical 
and numerical estimates show that a negative coupling constant is preferable for successful inflation. 
Moreover, provided fixed initial conditions, larger |^| leads to a greater number of e-folds. The latter 
is obtained for a reasonable range of initial conditions and the coupling parameter and indicates a 
possibility for successful inflation. 
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I. INTRODUCTION 

Modern astronomical observations gave rise to a num- 
ber of important conceptual issues, which can be success- 
fully resolved by the mechanism called inflation. Accord- 
ing to the inflationary paradigm 0, Q , there has to have 
been an epoch, when the Universe was undergoing an 
accelerated expansion; such that the size of the Universe 
increased by a factor of about e^" (60 e-foldings) to agree 
with observations. 

This scenario can be realized within the framework of 
the Friedmann-Robertson- Walker model, if the matter is 
represented by a self-interacting scalar field, the inflaton, 
which is minimally coupled to gravity. The mechanism 
considers so called slow-roll approximation, when the in- 
flaton rolls down its potential hill towards the potential 
minimum, but, at the same time, remains far away from 
the minimum. It is in the vicinity of the minimum of the 
potential, where inflation stops. For the simplest choice 
of a c^uadratic field potential with an arbitrary mass of 
the inflaton, one would obtain 60 e-foldings if the slow- 
roll regime starts off at « 3Mp and ends when the field 
is essentially zero. The question as to how the inflaton 
gets to that value before the slow-roll phase is left for a 
quantum theory to answer. 

Without a quantum theory taking into account gravity 
at hand, one could heuristically argue along the lines of 
Linde's chaotic inflationary scenario The idea was 
that vacuum fluctuations of a scalar field would be dis- 
tributed randomly up to Planckian energy densities. Re- 
gions with sufficient values of would inflate resulting in 
a flat FRW-uni verse. We shall see that, however insight- 
ful, this paradigm imposes severe restrictions on the like- 
lihood of inflation in non-minimally coupled models, with 
a coupling of the form — Specifically, the 
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coupling constant needs to be very small for successful 
inflation: ^ < 10^^^ for positive coupling and |^| < 10"'^ 
for negative coupling. 

Loop Quantum Cosmology (LQC) has over the 
last years resolved some long-standing issues of cosmol- 
ogy and general relativity in whole and suggested an al- 
ternative mechanism to chaotic inflation. As has been 
shown in 0, 0, LQC leads to effective quantum modi- 
fications of the classical Friedmann equations, such that 
gravity becomes repulsive at small scales and introduces 
an 'anti-friction' term into the Klein-Gordon equation 
for the scalar, thus driving the infiaton up its potential 

hiu i M- 

Therefore even small vacuum fiuctuations 
of (f> will be amplified during the quantum-corrected (ef- 
fective) phase of evolution. In the light of this, much 
smaller field initial values are required to start a success- 
ful classical inflation. This fact significantly relaxes the 
restrictions on the coupling constant of a non-minimally 
coupled theory. Moreover, it is possible to study the infla- 
ton's 'climbing-up' systematically, provided some initial 
fluctuations occur in the fleld and its canonically conju- 
gate momentum. 

The paper is organized as follows. We first review the 
conventional model for a non-minimally coupled scalar 
field. Then, following H Q, we perform a conformal 
transformation, which allows us to recast the action in 
the minimally coupled form in terms of a redefined field 
variable with some effective potential. In section three, 
we investigate the slow-roll regime and derive the number 
of e-folds as a function of the value of the inflaton right 
before inflation. Then, in section four, we turn to inves- 
tigating the inflaton's 'climbing', that is how small initial 
fluctuations get amplified during the quantum-corrected 
phase and the following classical evolution. The esti- 
mates for the maximum field value are mostly obtained 
analytically, up to a certain point when one is forced to 
solve a transcendental equation. The results are then 
compared with the numerical solutions of the Hamilton 
equations of motion in section five. The comparison indi- 
cates that the analytical method provides a decent preci- 
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sion and hence can be efficiently used for the analysis of 
the range of initial conditions and parameters that lead 
to a sufficiently long inflationary phase. 



II. ACTION FOR A NON-MINIMALLY 
COUPLED SCALAR AND CONFORMAL 
TRANSFORMATION 

The simplest form of a non-minimally coupled action 
is given by 



S[gab, 0] 



-9 (^R Ig'-'VacpV^cp 



(1) 

where g = dct(gab), n — SttG, R is the scalar curvature 
of the metric and U{(j)) is the self- interacting potential 
of the matter scalar field (j). The curvature term is cou- 
pled to matter through the coupling function f{(f>) and 
is manifestly Lorentz invariant. There are some obvious 
restrictions on the coupling function. First of all, as we 
shall see, it enters the symplectic structure for the gravi- 
tational (connection) variables. Thus it must not vanish 
anywhere. Moreover, one can view this function as defin- 
ing an effective gravitational constant k, in the sense 



cff = K/f{(p) 



(2) 



Hence we require f{(j)) to be always positive. Note also 
that the scalar field is expected to change sign, e.g. dur- 
ing the reheating epoch. It is then natural to restrict 
/((/)) to be an even function of (j). Finally, the coupling 
function must satisfy the limit f{4>) 1 when (jj —^ 0, 
so that the standard normalization is recovered in the 
absence of a scalar field. 

Given the restrictions above, the simplest form of the 
coupling function is 



1 — (7(7 



(3) 



(One can think of this expression as the Taylor expan- 
sion up to the quadratic order. We are interested in the 
case of a weak coupling, so the latter should be a good 
approximation. Note also that coupling of higher degree 
would require a dimensionful coupling constant, while ^ 
is unitless.) In fact, this is the form of coupling function 
that is normally considered in the literature (see, for ex- 
ample, 23). We will refer to a as the coupling strength 
(not to be confused with the coupling constant ^ = ^ of 
01 ) . Note that zero a corresponds to the case of minimal 
coupling. If (T > 0, the coupling function can, in prin- 
ciple, become zero and even negative if the value of the 
scalar field exceeds the critical value 



(4) 



We have deliberately written the absolute value of the 
coupling strength because, as we shall see later on, it is 
also of importance in the case of negative coupling. 



In this paper, we restrict our attention to a flat ho- 
mogeneous and isotropic model, described by the FRW- 
mctric 



ds^ = ^dt^ + a{tfdf^ 
Using a conformal transformation 
fiab ■■= .fi(f>)9ab 

the action can be recast in the form 



(5) 



(6) 



S[gab, ^] = dt 



(7) 



where we have introduced the effective potential 



y((/>) 



and 



1- (7(^^(1-^) 



(8) 



(9) 



The last approximation holds as long as we are interested 
in weak coupling such that a ~ (10"'^ 10~^)k; then 
(l — ^) can be set to 1 in this order of magnitude. Note 
that in the Hamiltonian formulation, there is a similar 
canonical transformation under which -F((/>)^ := ex- 
actly • Consequently, the relation between the scalar 
fields derived below, also becomes exact, which con- 
siderably simplifies the analysis even if ^ ^ 1. The last 
step, that will bring the kinetic term into its canonical 
form (Einstein frame), is to redefine the field 



d^F{<p) 



sinh"^( 



if (T > 
if (7 < 



(10) 



Note that for a positive a the quantity \/a(f) must be 
less than unity, for /((/>) must not vanish and, in fact, be 
positive. From H10|l we get 



am{y/a(f>), cr > 
sinh{y/\a\4>) , a < 



(11) 



Aiming at a loop quantization, we now have all the 
ingredients to proceed to the Hamiltonian formulation of 
the theory. For that end we introduce the phase space 
variables: 



■ 7^1 



and 



P 



(12) 



here p and c are the gauge-invariant triad and connec- 
tion components respectively, 7 is the Barbero-Immirzi 
parameter fl^ Il3 | , and tt is the field momentum. Gen- 
erally, the sign of p determines the orientation of the 
triad. From now on we assume that p > 0. Again, the 
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tildcd variables stand for quantities written in the Ein- 
stein frame. For notational convenience, however, we will 
omit tildes until transforming back to the Jordan frame. 
In terms of newly defined variables the action {Tj) takes 
the form 



dt 



— pc + 7r0 — Ti, 

K7 



(13) 



with the Hamiltonian 



n = — ^vpc' 



+ p"^V{4>) (14) 



From (|13|) we see that the pairs {p,c\ and {</>, tt} are 
indeed canonically conjugate variables with the Poisson 
brackets 



(15) 



Before proceeding to inflation, let us comment on a phe- 
nomenological relation between the Einstein and Jordan 
frames. More specifically, eventually we are interested in 
the number of e-folds during the inflationary epoch 



N ■= In 



at 



(16) 



where the subscripts stand for the initial (before infla- 
tion) and final (after inflation) values of the scale factor 
a. Suppose we have solved the equations of motion in 
the Einstein frame and obtained ah and Oe. Then it is 
clear from ^ that the scale factors in the two frames are 
related by 



a = VfWa 
Therefore the Jordanian number of e-folds is 



(17) 



= In 



V V/(0e)ab^ 



2 U(0e) 



(18) 

At the end of inflation, the scalar field is essentially zero, 
while at the beginning of inflation it has its maximum 
value (/)max- Thus, using we can explicitly write 



1 



N = -ln(l 



J 



(19) 



In principle, the beginning of inflation may occur when 
the scalar field is close to its critical value, such that 
Ic^maxl ~ 1. In this case, a negative coupling cr < 
would yield a negligible correction to the righthand side 
of H19() . A positive coupling, on the contrary, might give 
a significant contribution resulting in a substantial dif- 
ference between the Jordanian and Einsteinian number 
of e-folds. 



III. INFLATION 

In this section, we consider the inflationary mecha- 
nism 

M slow-roll approx^mat^onU■ 

According to this 

model, the inflaton rolls down its potential hill towards 
the potential minimum, but, at the same time, remains 
far away from the minimum. In such a regime, a scalar 
field behaves as a cosmological constant. It is in the 
vicinity of the minimum of the potential, where inflation 
stops. We compute the number of e-folds as a function of 
the initial (maximum) value of the scalar field. We shall 
see that the result is generic irrespective of the sign of 
the coupling strength a. 

We start with investigating the equations of motion 
generated by the transformed Hamiltonian (|14() . 



P = 




-j\lj + 2p^Vi^) 



(20) 



where the prime denotes a (/)-derivative. In the first 
equation we have eliminated the connection c, using the 
Hamiltonian constraint H14() « 0. We shall make an 
extensive use of these equations later on. For the mo- 
ment, however, it is more convenient to consider the con- 
ventional Friedmann equations. The constraint equation, 
rewritten in terms of the scale factor and its time deriva- 
tive, reads 



^„^2 ^ « ^2 ^ ^3-^(^) ^ Q (21) 
K 2 



n = 

Dividing by a^, one gets the Friedmann equation 



2 



(22) 



Here H — a/ a is the Hubble rate. Using the 7r-equation 
we derive the Klein-Gordon equation 



(23) 



Finally, the Raychaudhuri equation can be obtained by 
taking the time derivative of l|22f) and substituting </> from 
(EH): 



H 



(24) 



From now on we stick with the slow-roll approximation 
that would ensure inflation. The potential domination 
(small kinetic terms) 



</.2 < 1/(0), \4>\ < \3H^\ 



(25) 
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is provided, if we require the slow-roll parameters to be 
small; 



(26) 



With the assumptions (|25|) the Friedmann and Klein- 
Gordon equations can be rewritten as 



3i7^« -¥'{(/)) 



(27) 



V(x) 




whereas the Raychaudhuri equation simply implies H w 
0, i.e. an exponential expansion of the scale factor. Di- 
viding the first of the equations (|27|l by the second one, 
we can express the Hubble rate as 



H : 



(28) 



The number of e-folds is then given by 



N 



Hdt: 



V 



(29) 



where ^f, and are the values of the inflaton before and 
after inflation respectively. The last expression tells us 
that as soon as the shape of the potential is specified, one 
can compute the number of e-folds. Therefore, it is now 
pertinent to analyze the form of the effective potential 
given by (jHJ. The most common choice is either quadratic 
or quartic potential 



-02, or C/(0) 



A 



2 " ' v^/ 4 

The corresponding effective potentials are 



(30) 



F(0) = 
if a is positive, and 

y(0) 



2(7 cos^(-ya-0) ' 
A 
4^ 



and 



tan^(v^'/') 



(31) 



tanh^(v/|cr|0) 
2|(t| cosh2(yH0) 



and 



y(0) = ^tanh4( 



(32) 



if a is negative. 

Fig. n shows the effective potentials for the quadratic 
original potential. If a > (Fig. the effective poten- 
tial is steeper than quadratic (punched line) and diverges 



FIG. 1: The effective potential for a < (upper curve) and 
(7 > (lower curve). The solid curve represents the original 
quadratic potential. The argument of the potential is x = 
\f\o\(i) for (T 7^ and x = (f) for cr = 0. 



at \/a4i = ^ thus keeping the inflaton away from that 
point. Note that this is exactly the value of at which 
the coupling function goes to zero. On the contrary, the 
effective potential for a negative a is bounded from above, 
goes below the quadratic potential (punched line) and has 
a maximum at -y/joj^ = ln(l + ^/2). Remarkably, both 
these values of the tilded scalar correspond to the crit- 
ical value of the original field Q). The existence of the 
maximum for a negative coupling constant implies that 
if the inflaton exceeds the critical value, it will roll down 
to the right, towards infinite values. The latter is unac- 
ceptable, as the inflaton is supposed to dissipate to zero 
after inflation. We therefore should restrict our attention 
to < 0crit for both positive and negative couplings. 

Let us now retrieve the tildes over the conformally 
transformed variables, to distinguish them from the orig- 
inal ones, and calculate the number of e-folds 



N 



V 

— — ( 

01. ^.^ 



(33) 



1 



d4) 



0. /(</*) In fmV 
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If the original potential is quadratic U {(f) ~ ^m'^(j>'^, then 



(123 yields 



N : 



1 + cr<7 



4(7 



■In 



1 +(70 



(34) 



Neglecting (j)^ as compared to (/>{, = (t>ms^x, and using Hll|l . 
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the number of e-folds in the Einstein frame is 



A. Effective phase 



N w — In 1 

4(7 ' 



K 

4^ 



ln|l + sin^(\/CT0niax^ 

hi|l -sinh2(yH0n 



cr > 
cr < 



(35) 



Note that in the hmit ct ^ 0, the number of e-folds re- 
duces to that of the minimaUy coupled scalar field. Fur- 
thermore, considering the Taylor expansion of H35|) 



-ln|l 

(7 



2^4 

X ax 

2 



and the inequality | sinh(a;)| > | sin(a;)|, we conclude that, 
provided the same initial field, negative cr is preferable for 
successful inflation. 

Following a similar derivation for a quartic original po- 
tential U{(p) = |-</>^, we obtain the number of e-folds 



N : 



sin^(VCT^niax), 

sinh^(y|CT|^ma, 



cr > 
(T < 



(36) 



As before, in the limit cr ^ the expression agrees with 
that of the minimally coupled model. Furthermore, a 
negative coupling constant again provides a greater num- 
ber of e-folds. 

As was mentioned before, the question as to how the 
inflaton gets at 0max is not addressed in the classical 
theory. In the next section, we shall estimate (/)max on 
the grounds of LQC. 



IV. ESTIMATION OF MAXIMUM INFLATON 
VALUE 



In the framework of LQC, the self interacting scalar 
field, infiaton, arises as a microscopic vacuum fluctuation 
that is driven up its potential well by quantum effects. 
In the isotropic and homogeneous context, it is indeed 
quantum corrected (effective) Friedmann equations 0, 1^ 
that have been shown to govern the evolution of the Uni- 
verse in the semiclassical regime, when the scale factor 
ao < a < a* . The boundedness of the geometrical density 
operator manifests itself in the form of "repulsive" 
effects of gravity at Planckian scales and leads to an 
anti-friction term in the effective Klein-Gordon equation 
i 

Practically, the process of the inflaton's 'climbing- 
up' can be broken into two stages: i) effective, driven 
by quantum modifications and consecutive ii) classical 
phase, conditioned by the field momentum gained in 
the effective regime. Both stages will be assumed to be 
kinetic- dominated, that is they occur relatively rapidly, 
and the potential terms, whenever they appear along 
with kinetic ones, are dominated by the latter. In other 
words, (j) almost reaches its maximum before the poten- 
tial effects become appreciable. 



This phase starts with some initial conditions and ends 
when the scale factor a ~ a, , a characteristic scale intro- 
duced in Ea. (|39|l below. We should say that the idea 
behind the estimate for TV used in this paper is analo- 
gous to that in except that the semiclassical phase 
ends before the inflaton reaches its maximum. Nor do we 
assume any specific asymptotic limit of the geometrical 
density spectrum. 

The recipe provided by LQC as to how to proceed 
directly to an effective description is the following. In 
the Hamiltonian constraint, one is to replace all negative 
powers of triad p with appropriate factors of the spec- 
trum of the inverse volume operator [T^ 



where 



Dig) = q 



— 



1 

^372 
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d,= 



2,3/2 



(37) 



(q + l)"/^-|g-l 



11/4 



(<Z + l)^/4-sgn(g-l)|g-ir/4]}'(38) 



and 



_ P 2 87r7j/io 2 ^on^ 

q = — with = := Ip (39) 

P* 3 



The Planck length is defined as £p = ■^/kTStt = AIp^. 
The asymptotic behavior of the density operator fol- 
lows from the equation above. Not only does dj remain 
bounded for p ^ 0, it becomes proportional to a positive 
power of the scale factor. Specifically, D{p/p^) oc p^'^/^ 
and d (x p^, thus making the matter Hamiltonian H38|) 
well behaved. On the other hand, as p goes to infin- 
ity, one recovers the classical dependence: D{p/p^:) w 1 
whereas d ~ p~^/'^. 

Note, however, that this Hamiltonian constraint does 
not account for higher order perturbative corrections that 
proved to be important for small values of the parameter 
i Hi ^1 More precisely, one has to consider the 

modified Hamiltonian 

Hcff = ^VPsin2(^oc) + ^dj(a)^' + V{(j,) (40) 

When hqc <C 7r/2 one recovers the Hamiltonian H14|l . 
The discreteness corrections become important, when the 
matter part of the Hamiltonian is of order of the critical 
density 



o / 222 
Per = 3/kpo7 a 



(41) 



corresponding to the maximum value of the gravitational 
part, attained when sin(/xoc) = 1. 

Let us focus on the case, when the initial conditions are 
such that /xpc ^ 7r/2. We shall investigate under what 



6 



restrictions this inequality will hold during the evolution 
within the effective regime. The effective counterparts of 
equations of motion H20(l take the form 




D{q) — + 2pW{(t)) 



^ ^ D{q) 



,3/2 



(42) 



As usual we have eliminated the connection variable us- 
ing the constraint equation 



c 

7 



(43) 



Note that within the effective regime, the field is very 
small {a4>^ 1) the non-triviality of the coupling 
can be neglected. Let us therefore set /(0) to unity. 
From now on we take the quadratic potential 



1 



2 J,2 



V{(j)) = -TO>- 



(44) 



m being the field mass. Consider the initial stage of the 
effective evolution, when the scale factor ranges from p w 
Po = ^^^3^° to p = p*. For a small value of the inflaton 
((f) being of order 0.1), we see that in (|43|l the kinetic 
term overwhelmingly dominates over the potential term, 
provided the field mass m = 10~^Mp. Thus for the sake 
of estimate of the maximum of ^qc we can discard the last 
term. Moreover, we notice that multipliers in the first 
term can be split into "fast" and "slow" ones. Formally, 
TT is a slow variable, as the righthand side of Ti-equation 
in (|42|l is small due to a factor of m? in the potential 
term, whereas p is a fast variable. As far as D{q) is a 
large positive power of p, the D{q)-i&ctoT in H43|) is the 
fastest. More rigorously, this statement can be cast as 



D{q) 



D{q) 



> 



> 



(45) 



Although, (j) may become as fast as p, especially in the 
end of the effective phase, it is not used in the following 
derivation. The time derivative of c can be approximately 
written as 



7 V 6dt Vp 



This expression vanishes when 



K7rd(/D(^) 
6 p 



dt 



(46) 



dp 



which is equivalent to 



dD{q) 
dq 



The solu tion to t he last equation is (?max ~ 0.97, which 
yields V^(9max)/gmax « 1-2. Neglecting the change in 
TT and taking its initial value, we obtain the estimate for 
the maximum value of poc 



K \/D{qjnax) TTi 0.3 

^ fJ-Oj— = — X TTj (47) 



This remarkable result tells us that, provided the afore- 
mentioned restrictions, the maximum value of /iqc, and 
hence the maximum matter density, depends only on the 
initial field momentum, regardless the initial scale factor 
and field itself. Specifically pmax is linearly proportional 
to TT?. Thus the approximation ^qc <C 7r/2 remains valid 
throughout the effective regime, if 

^^^ < 5j£p (48) 

Note that a similar derivation goes through, if the modi- 
fied Hamiltonian (|40|) is used. At a given scale factor, the 
matter density must not exceed the critical value H41(l . 
The condition on the initial momentum H48|) then be- 
comes 



TTi < 



2k 



3D{q„ 



3.45j4 



This condition is not very restrictive and compatible with 
characteristic values of initial field momentum fluctua- 
tions derived from the Heiscnberg uncertainty principle, 
provided initial held fluctuations of order 10^^i\/p. For 
a more detailed analysis see [T^ . 

Now suppose that the effective evolution starts at some 
initial data p = pi,(f) = (j)i,TT = TTi. If 'I'i satisfies the con- 
dition H48|) . we can safely use the approximate Hamilto- 
nian (|37|) and Hamilton equations H42|l . In order to esti- 
mate the maximum value of the inflaton at the end of the 
effective phase, it is convenient to rewrite the equations 
above in terms of the following dimensionless variables 



{mt)^t, X := — , ?/:=—, 
P* Alp 



mMppl^'^ 



(49) 

In these variables the equations of motion 142(1 take the 
form 



IQtt L^, .z'^ t , 
3 V X 
= D{x)- " 



.3/2 



(50) 
(51) 
(52) 



We are going to get the approximate values for y and 
z (hence and tt) at the end of the effective phase. In 
the same fashion, as was done before, we can integrate 
Eq. (|51|l . treating z as a slow variable and fixing it at a 
constant value, that corresponds to Zi. Dividing H51|l by 
(I50|l . we express the derivative 



dx 



D{x) 



(53) 



7 



then the integration yields 



y{x) ~yi + \ hr / ; ax . 

V 2k L.. x' 



(54) 



the right-hand side of the 7r-equation. In the hght of this, 
we eliminate the time derivatives in Eqs. H2U|I : 

dp 



Similarly 



dx 



3 x'^y{x) 



(55) 



Where we again have neglected the potential term in (|5 
Integrating we get 



z\x) 



Direct calculation shows that given reasonable initial con- 
ditions, the changes in the field and its momentum are 
indeed small and can be neglected when computing the 
value of the maximum matter density. In the dimension- 
ful variables, the value of the inflaton at the end of the 
effective regime, 0Oj is 



dq 



(57) 



Assuming qi = —j 
scalar field is 



0.1 and go ~ the change in the 



0.14Mp 



(58) 



The subscript 'eff' stands for effective. Remarkably the 
growth of the inflaton does not depend on its initial mo- 
mentum at all! 

We should now take the obtained values for po and 00, 
as well as ttq ~ tt^, to be the initial data for the classical 
phase and estimate how high the inflaton will get. 



B. Classical phase 

When the scale factor significantly exceeds a*, the 
quantum corrections become negligible and the evolu- 
tion is very well described by the classical equations (|20|l . 
Note, however, that we can no longer consider the infla- 
ton to be small and should use the effective potential. 

The end of the 'climbing' phase, i.e. the beginning of 
inflation, occurs at cx tt = 0, which docs not corre- 
spond to any a-priori fixed scale factor. This is exactly 
the reason why one cannot use the technique of the pre- 
vious sub-section. Instead one ca n ag ain assume kinetic 
domination almost up to the top [IJ, where (j) = (/)max- 
In other words, the inflaton changes insignificantly af- 
ter the kinetic terms become less than the potential ones 
near </)max- Note that we neglect the potential only when 
it appears along with a kinetic term; we should still keep 



do TT 



(59) 



Solving the first equation for the scale factor as a function 
oi(b 




P{4') = Po cxp { M 



(60) 




and substituting this expression into the second equation, 
we obtain 



dn 



^exp{V6^(0-0o)}l^'(0) (61) 



Finally, separating the variables in the last equation and 
integrating it from tt = ttq (beginning of the classical 
climbing) to tt = (end of the classical climbing), we 
arrive at a transcendental equation for (/)max: 



2pg 



exp 



4>o 



{V6^i(j)-M}v'Wd(b (62) 



Before giving numerical solutions, let us draw some an- 
alytical conclusions from this equation. As the effective 
potential is steeper for a positive a, we should expect a 
smaller maximum value of the inflaton in this case than 
for a negative coupling. Note also that in the integrand 
of H62|) the exponential is the fastest growing function, in 
the same sense as before. Therefore, when integrating, 
the potential factor can be taken out of the integral and 
evaluated, for the sake of estimate, at the upper limit. 
We now restrict our attention to the case of quadratic 
original potential and consider the Taylor expansion of 
the effective potential 

F(0) = imV(l + ^^'^' + O(04)) (63) 

Then the transcendental equation (|62|l can be approxi- 
mately rewritten as 



a 

K 



where we have defined 
y := V6K{(f)m 



A 



9 T 



(64) 



(65) 



Let us first set a to zero and denote the corresponding 
solution y, i.e. y satisfies 



ye^ = A 



(66) 



The solution to this equation is the Lambert function 
y = Ty(A). It is plotted in Fig. |21as a function of the 
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1.8 ^ ' ' ' 1 
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JTof/p) 

FIG. 2: Maximum change of the inflaton as a function of the 
initial field momentum for cr = and po = p* ■ The analytical 
curve is the Lambert function defined in 166L whereas the 
numerical points are obtained as solutions to the equations of 
motion lH^ . 

initial field momentum for a fixed initial scale factor po ~ 
p,. We see that dependence of 0max — 0o upon ttq is 
close to logarithmic, i.e it is not very sensitive to initial 
momentum. For instance, if tto ~ 5£p, the inflaton grows 
by about two Planck masses during the effective phase. 
This is insufficient for successful inflation. 

Let us now take a non-zero coupling constant and, 
treating it as a small parameter, consider the corrections 
to the maximum value of the scalar field e := y — y. The 
equation H64|) then yields 

<^ = -7ry ^ 5-CT(0max - (po) (67) 

The first important consequence of this expression is that 
the fT-corrections to ^max are small and proportional to 
minus a. In other words, in the case of a negative cou- 
pling strength, the inflaton indeed gets higher up the 
potential hill. 

Before combining all the contributions to (/>max, we note 
that there also takes place an effective (super-) inflation, 
while the universe is expanding from a — ai to a ^ 
0- This also provides a couple of e-folds, as the scale 

factor a = y/p increases by a factor of = ■ 

According to ||2SJ), (pmax ~ 2.6Mp is sufficient to pro- 
vide 60 e-folds. Let us now give a numerical estimate for 
the maximum inflaton value. For j = 5, taking the ini- 
tial conditions Pi = p*/2j = O.lp*, 4>i = O.lMp, tt^ = 5£p, 
the effective growth of the field is O.I4A'/p, whereas the 
classical contribution is 2.25Afp. Adding everything up 
we get 0max ~ 2.5A'/p. We shall discuss how generic this 
result is and how it varies with a different choice of initial 
conditions in the next section. 



V. NUMERICAL RESULTS AND DISCUSSION 

In this section, we study numerical solutions to the 
differential equations (|42|) . Again we will be mostly in- 



terested in the maximum value of the inflaton for the 
allowed range of initial conditions and the parameter a. 

Let us first discuss the initial field and the scale factor. 
As was shown in the previous section, the dependence on 
the former is rather trivial: the initial value 0o acts as 
constant of integration and should be merely added to 
the change in the inflaton during the climbing phase. 

The dependence of 0max on the initial scale factor is 
more complicated, but the value of pi itself is quite re- 
stricted. It appears that there is a natural choice of pi, 
associated with the smallest eigenvalue of the area oper- 
ator [23, related to /xq. In fact, as one can see from ((STJ, 
the effective growth of the scalar A0off depends on the 
limits of integration, hence scale factor, only through the 
ratio q=^. The lower limit is given by g,; = while 
the upper limit is not fixed a-priory. It has a physical 
meaning of the marginal value that separates the effec- 
tive and classical behavior of the spectrum of the geomet- 
rical density operator. For the estimate in the previous 
section we used g = 1 as the upper limit. A closer look 
at the graph of D{q) (see Fig. O shows that there is a 
maximum around q k. \. In other words, the behavior 
of the spectrum is not yet classical. It would be more 
reasonable to stop the integration (which is meant to be 
over the effective domain) at a somewhat greater value 
of (7, where D((i) is essentially close to unity. Recall that 
go = 1 was giving A(/)off = 0.14Mp. At the same time, 
analysis of Eq. H54() , for A^eff as a function of the upper 
limit of integration goj shows that the effective growth of 
the inflaton sensitively depends on q^, whereas the clas- 
sical formula Ht)t)|l indicates a much weaker dependence 
upon go- For instance, taking go = 2, one would get 
A(/)cff = 0.3A/p, while go = 4 yields A^off = O.SAfp. 
This implies that if one increases go, one should expect 
a somewhat greater value of total A(/) = A^cft + ^4>c\- 
Nevertheless, for the sake of estimate, we will stick with 
go = 1 and compare the numerical results with analytical 
formulas. 

Summing up, as the effective growth of the inflaton is 
not very substantial and does not depend on the initial 
field momentum tto or coupling strength cr, the quantity 
of the most interest is the classical part of A0. The latter 
is determined by tto and a. 

The dependence of A(/)ci on the initial momentum 
Po for five different values if a is displayed in Fig. 
^ The curves correspond (top to bottom) to cr/-^p = 
—0.1,-0.05,0,0.05 and 0.1 and are plotted for the al- 
lowed range of initial momenta: |^ < 5j' = 25. The 
maximum inflaton value grows with initial momentum 
and qualitatively resembles the Lambert function. It is 
steepest for small momenta and flattens out when ttq be- 
comes large. Quantitatively, for cr = there is a very 
good agreement with the analytical result H66(l and the 
graph of Fig. [3 The curves, corresponding to opposite 
values of the coupling strength, are situated symmetri- 
cally around the minimal curve, which justifles the linear 
(in cr) approximation used in H63|) . Furthermore, direct 
calculation shows that the cr-corrections to the central 
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FIG. 3; The spectrum of the geometrical density operator 
D{q) introduced in Eq. (1^ 
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FIG. 4: The change of the inflaton during the classical phase 
as function of its initial momentum for different values of a 
(measured in the units of £p). 

curve are very close to those predicted by formula H67() 
for both small and large initial momenta. 

Classical climbing brings the scalar field fairly close 
to the required values determined by iV = 60 and (|35|l . 
For instance, (p^cq = 2.6Mp for cr = — 0.05^| and (/freq = 
2.5Mp for a = — 0.10.^|,). However, A0ci on its own is 
not enough for sufficient inflation. At the same time, the 
number of e-folds is very sensitive to deviations in 0max 
near the required values. In other words, even small (pos- 
itive) variations of </>,; and/or A0cff, such as of order of 
O.lAfp, would lead to sufficient and should be seriously 
taken into account. 

We see from Fig. 0] that a negative coupling con- 
stant indeed gives a greater change of the inflaton, and 
more negative values are preferable for successful infla- 
tion. This fact can be also understood from the Jor- 



danian perspective. The non-minimally coupled action 
(Q) can be thought of as a standard one with a field- 
dependent gravitational constant ((SJ. Clearly, a nega- 
tive cr results in 'weaker' gravitational coupling, and it is 
not surprising that the inflaton would reach to a greater 
value (yielding a greater N), than for <t = 0. Moreover, 
if (7 < 0, the second term of the righthand side of the re- 
lation between the number of c-folds in the Einstein and 
Jordan frames (|19|l is of order one and can be neglected. 
Thus N N and is greater than in the case of mini- 
mal coupling. Similar considerations work for a positive 
coupling as well and imply a smaller number of e-folds. 

We should now clarify the seeming discrepancy with 
the result of Futamasc and Maeda 0| . As wc have already 
mentioned, in that paper, the authors argued that, in 
order to allow successful inflation, the coupling constant 
^ had to be fine-tuned within a very narrow interval: 
^ < 10-^2 for ^ > and |^| < lO^^ for ^ < 0. These 
restrictions on ^ arose from the heuristic argument based 
on Linde's chaotic inflationary scenario. The values of (pi 
arc assumed to be randomly distributed from zero up to 
Planckian energy densities V{(f)) ~ Mp. Such potential 
is attained at ~ lO^Afp. The latter must necessarily 
be less than the critical value cpd-it Q, which implies 
a ~ 10-124 and C ~ 10"". 

The above conditions can be relaxed in the framework 
of LQC. As the initial values of the scalar field appear 
as vacuum fluctuations, that are amplified during the 
'climbing' phase, one just needs to restrict the coupling 
constant so that 0crit merely exceeds (j) ^ 3Mp - the max- 
imum inflaton value we are interested in. That yields 
a ^ 10~^£p and ^ ~ lO""^. The former is exactly the 
maximum value of the parameter a we have considered 
in the paper. 

To summarize, the main characteristic of inflation, the 
number of e-folds, depends on the initial matter fluc- 
tuations and the coupling parameter a. The coupling 
strength is bounded to be less than 0.1£p by the condi- 
tion (/)max < (Acrit = 1 / 'v/H' both positivc and neg- 
ative couplings. At the same time, the most negative a 
would work best for successful inflation. The restriction 
on the initial field momentum tt^ appears as a require- 
ment for the matter density to remain subcritical and 
implies |7r,i| < 5j. With these restrictions, the value of 
the inflaton will increase by approximately 2.0 — 2.6Mp 
during the 'climbing' (effective and classical) phase. To- 
gether with the initial vacuum fluctuations of the scalar 
field of order of several tenths of Mp this would lead to 
N > 60, i.e. sufficient inflation. It should be noted that 
one cannot get N much greater than 60, which indicates 
that observations may be sensitive to the mechanism dis- 
cussed here. 
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